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Gravitational-radiation losses from the pulsar— white-dwarf binary PSR J1141— 6545 
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Pulsars in close binary orbit around another neutron star or a massive white dwarf make ideal 
laboratories for testing the predictions of gravitational radiation and self-gravitational effects. We 
report new timing measurements of the pulsar-white-dwarf binary PSR J1141— 6545, providing 
strong evidence that such asymmetric systems have gravitational wave losses that are consistent 
with general relativity. The orbit is found to be decaying at a rate of 1.04 ± 0.06 times the general 
relativistic prediction and the Shapiro delay is consistent with the orbital inclination angle derived 
from scintillation measurements. The system provides a unique test-bed for tensor-scalar theories 
of gravity; our current measurements place stringent constraints in the theory space, with a limit 
of Qq < 2.1 x 10~ 5 for weakly non-linear coupling and an asymptotic limit of Qq < 3-4 X 10 -6 for 
strongly non-linear coupling, where «o is the linear coupling strength of matter to an underlying 
scalar field. This asymptotic limit is nearly three times smaller than the Cassini bound (qq ~ 10 -5 ). 

PACS numbers: 04.30.Tv, 04.80.Cc, 95.30.Sf, 95.85. Sz 



Introduction. — Einstein's general theory of relativity 
(GR) has passed all experimental tests so far with com- 
plete success, which makes it one of the most celebrated 
theories in modern physics [l|. While the most precise 
tests have been conducted in the weak-field conditions 
of the solar system @, y], massive and compact astro- 
nomical objects such as neutron stars and black holes, 
in particular pulsars in close binary orbits with another 
neutron star or a massive white dwarf (i.e. relativistic bi- 
nary pulsars), allow testing GR in strong-field conditions 
0; El @ • These tests have confirmed GR at an impressive 
level of better than 1% @, i, i, i, @ . 

While GR is indeed the most successful theory of grav- 
ity, it may conceivably break down under extreme strong- 
field conditions where other theories of gravity may ap- 
ply. Tensor-scalar theories, which invoke a coupling be- 
tween matter and an underlying scalar field (in addi- 
tion to the standard space-time tensor field), are thought 
to be the most natural alternatives to GR [7], Q . The 
most well-known in this framework is the Jordan-Brans- 
Dicke theory, with a linear coupling between matter and 
scalar field: a(ip) = ao ip, whilst a more general de- 
scription involves a two-dimensional parameter space: 
a(tp) — ao"0 + \flo V' 2 > where tp is the strength of the 
scalar field; a(ip) is the coupling strength between the 
scalar field and matter; ao and Po are the coupling param- 
eters. The scalar field does not exist in GR and therefore 
(ao , (3o ) = (0, 0). Possible deviations from GR are thus 
constrained by experimentally imposed bounds on ao and 
3n . While several tests have been devised for this (see 
[a IS LLQ| for recent reviews), current best limits come 
from timing binary pulsars (/3o > —4.5) and the Cassini 
time-delay experiment (a^ < 10~ 5 ; Q). 

The binary-pulsar tests so far d, [5J, Q have focused on 
systems that consist of two almost identical neutron stars 
(i.e. symmetric systems in which both the objects are of 
nearly equal mass and size). PSR J1141— 6545, on the 
other hand, comprising a strongly self-gravitating neu- 
tron star and a (relatively) weakly self-gravitating white 



dwarf, is a gravitationally asymmetric system and thus 
provides a unique laboratory for testing GR and alterna- 
tive theories of gravity [8| . The asymmetry here is due to 
the self-gravity (e) or the compactness of the two bodies 
(given by e = —GM/(Rc 2 ), where M and R are mass 
and radius respectively). Since e ~ —0.2 for a neutron 
star and ~ — 1CP 4 for a white dwarf, a neutron star- white 
dwarf binary is a very asymmetric system. 

Discovered in 1999, PSR J1141-6545 is a relatively 
young pulsar (characteristic age ~1.4 Myr) spinning at a 
rate of once every 394 ms and is in a 4.74-hr orbit with 
a moderate eccentricity of ~0.17 [ll|. Early studies sug- 
gested that the pulsar lies at a minimum distance of 3.7 
kpc [12] and the orbit is inclined at an angle of 76° ± 2°. 5 
with respect to the plane of the sky [13]. The orbital 
period derivative from initial timing analysis was shown 
to be consistent with GR at the 25% level [14] and, more 
recently, observations over a 6-yr time span have unveiled 
remarkable changes in both the pulse shape and polarisa- 
tion, which are attributed to geodetic precession resulting 
from general relativistic spin-orbit coupling 15]. 

In this article, we report measurements of four post- 
Keplerian parameters for PSR J1141— 6545: advance of 
periastron to , time dilation and gravitational redshift 7, 
the Shapiro delay parameter s, and the orbital period 
derivative P\, , that allow us to determine the masses of 
the pulsar and its companion; demonstrate that gravi- 
tational wave radiation losses are consistent with those 
predicted by GR; and place stringent limits on tensor- 
scalar theories of gravity. 

Timing Observations and Analysis. — Timing observa- 
tions of PSR J1141-6545 were undertaken using the 64- 
m Parkes radio telescope in New South Wales, Australia, 
between 2001 January and 2007 May. Data were recorded 
primarily in the 20 cm radio band, while limited data 
were gathered in the 50 cm band. Early observations 
(2001 to 2002) made use of a multichannel filterbank 
as the backend, recording data over a bandwidth of 256 
MHz, while data from 2003 onward were recorded us- 



TABLE I: Timing parameters of PSR J1141-6545. 

Spin and astro-metric parameters: 

Right ascension, a (J2000) ll h 41 m 07 s .0140(2) 

Declination, 8 (J2000) -65°45'19". 1131(15) 

Pulse frequency, v (s _1 ) 2.538729369926(1) 

Reference epoch (MJD) 51369.852500 

Dispersion measure (pccm -3 ) 116.080(1) 

First derivative of pulse frequency (s~ 2 ) -2.767986(1) xl0~ 14 

Keplerian parameters: 

Orbital period, P h (days) 0.1976509593(1) 

Projected semi-major axis, x (s) 1.858922(6) 

Orbital eccentricity, e 0.171884(2) 

Time of periastron passage, T (MJD) . 51369.8545515(9) 

Longitude of periastron, u> (°) 42.4561(16) 

Post-Keplerian parameters: 

Advance of periastron, Co (° yr" 1 ) 5.3096(4) 

Gravitational redshift parameter, 7 (ms) 0.000773(11) 

Orbital period derivative, A -0.403(25) x 10~ 12 

Figures in parentheses represent the nominal 1-cr 
uncertainties in the least-significant digits quoted. 



ing the Caltech-Parkes-Swinburne baseband recorder 2, 
which coherently dedisperses data over two adjacent 64- 
MHz bands, thus yielding a net bandwidth of 128 MHz. 
The data were folded at the predicted topocentric period 
of the pulsar, where the adopted integration time var- 
ied from 60 seconds (2001 to 2003) to 300 seconds (2004 
onward), to account for almost a factor of two decrease 
observed in the pulsar's mean flux density over this pe- 
riod. The measurement uncertainties in the pulse arrival 
times have significantly degraded over this period, as a 
result of a dramatic increase observed in the pulse width 
(nearly by a factor of two over seven years) probably 
caused by precession of the pulsar spin axis. 

The pulse arrival times were computed by correlating 
the observed pulse profiles with template profiles of high 
signal-to-noise ratio constructed from long integrations 
of data. A total of 12,842 pulse times-of-arrival (TOAs) 
were measured. The remarkable changes observed in the 
pulse profile over our seven-year long observing span, 
while exciting for studies of geodetic precession and mod- 
elling the pulsar emission geometry, have complicated our 
timing analysis. In order to minimise the systematics 
caused by secular profile changes, we adopted a strategy 
which involves the use of a standard (template) profile 
that is a function of time. The final TOAs were anal- 
ysed using the standard timing package TEMPO, fitting for 
the pulsar spin, astrometric and Keplerian parameters, 
as well as three post-Keplerian parameters according to 
the relativistic and theory-independent timing model of 
Damour and Deruelle (DD) [16}. The final root mean 
square post-fit residual is 154 /is. The measured pulsar 
and binary system parameters are listed in Table 1. 

Tests of General Relativity. — In double neutron star 
systems and in close eccentric systems such as PSR 



J1141— 6545, the gravitational fields are so strong that 
the application of relativistic gravity becomes essential 
in timing models. For such systems, the observed pulse 
arrival times are modified by relativistic effects, which are 
potentially measurable through long-term (several year) 
timing observations. These relativistic effects may man- 
ifest themselves in various ways; for instance a temporal 
change in the period or orientation of the orbit, or an ad- 
ditional time delay (the Shapiro delay) resulting from the 
curvature of space-time when pulses pass near the mas- 
sive companion. These effects can be modelled in terms 
of the post-Keplerian (PK) parameters, which are essen- 
tially some phenomenological corrections and additions 
to the Keplerian orbital parameters [Tg|. These PK pa- 
rameters have different dependencies in different theories 
of gravity, and thus facilitate important tests of theories 
of gravity [l], [l7| ■ In any theory of gravity, the PK pa- 
rameters can be expressed as functions of the pulsar and 
companion masses and the easily measurable Keplerian 
parameters. A binary pulsar requires two PK parame- 
ters to completely determine the inclination and masses 
of the system. Measurement of three or more PK param- 
eters over-determines the system, and thus can provide 
one (or more) test(s) of gravitational theories through 
self-consistency checks. 

The orbital period derivative Ph is a crucial parameter 
as it is still the only observable that has ever verified 
the existence of gravitational waves. Relative to pre- 
vious timing analysis [14j . the measurement precisions 
on the PK parameters have now improved by up to a 
factor of 4 and, most notably, Ph is determined to be 
(-4.03 ±0.25) x 10~ 13 . This measured value (P^ cas ) is 
a combination of the orbital decay due to the emission 
of gravitational radiation (P^ R ) and contributions re- 
sulting from both real and apparent accelerations of the 
binary system along the line of sight; i.e., 



pmcas pGR 



pkin 



pGal 



(1) 



where "Gal" and "kin" refer to the Galactic and kine- 
matic contributions respectively [181 ] . The kinematic con- 
tribution, known as the Shklovskii effect (an apparent 
acceleration of the system due to its space motion, v t ), 
is given by vf/(cd), where d is the pulsar distance and c 
is the speed of light. The Galactic contributions include 
differential rotation in the plane of the Galaxy and accel- 
eration in the Galactic gravitational potential. Given a 
transverse speed of 115 ± 15 km s _1 deduced from scintil- 
lation observations and a pulsar distance of 3.7 kpc, we 
estimate P b kin = (6. 71±1. 73) xl0~ 15 and P b Gal = (-5.05± 
0.44) x 10~ 15 . Subtracting these from the measured P\, , 
we obtain an intrinsic value (—4.01 ± 0.25) x 10 -13 for 
pGR _ xhis is 1.04 ± 0.06 times the general relativistic 
prediction and corresponds to a shrinkage of the pulsar's 
orbit at a rate of approximately 2 mm per day. 

As the measurement of any two PK parameters allows 
solving for the two unknown stellar masses, our mea- 
surement of three PK parameters offers an independent 
test of GR. These results can be displayed elegantly in 



a "mass- mass" diagram as shown in Figure 1. Measure- 
ment of the PK parameters gives curves on this diagram 
that are, in general, different for different theories of grav- 
ity but should intersect at a single point (i.e. at a pair of 
mass values) if the theory is valid. Our results confirm 
that GR is a correct theory of gravity for asymmetric bi- 
nary systems, thereby extending the range of systems for 
which GR provides an accurate description. 

Given this, it is justifiable to apply the Damour Deru- 
elle GR (DDGR) formalism (as implemented in TEMPO) 
to determine the pulsar and companion masses [la . Il7j ] . 
This timing model is in the framework of DD but assumes 
GR to be the true theory of gravity. It uses measurements 
of the PK parameters u and 7 to solve for the companion 
and total masses of the system (m c and m to t respectively) 
and to model the Shapiro delay. The pulsar mass (m p ) is 
then given by m to t — m c • The best-fit values from such a 
DDGR fit to our timing data are m tot = 2.2892 ± 0.0003 
M Q and m c = 1.02 ± 0.01 M , which implies the pulsar 
mass TO p = 1.27 ± 0.01 M Q . These mass estimates im- 
ply an orbital inclination angle (i) of 73° ± 1°.5. This 
is in excellent agreement with independent constraints of 
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FIG. 1: Graphical summary of relativistic parameters from 
timing measurements of the PSR J1141— 6545 binary sys- 
tem. The shaded area is excluded from a consideration of 
Kepler's laws and the other constraints are based on the mea- 
sured PK parameters (shown as pairs of lines, with the line 
separation indicating the measurement uncertainty) inter- 
preted within the framework of GR. The systemic mass mtot = 
2.2892 ± 0.0003 Mq is very accurately defined from Co and 
is shown by the straight solid lines; the 7 term constrains 
the system to lie between the dashed lines, and P a between 
the pair of curved solid lines. The sin % measurement from 
scintillation analysis [13| is shown as the pair of dot-dashed 
lines. The hatched areas (inset) accentuate the masses al- 
lowed within the framework of GR. 



80 


I 


1 - 


60 






: 












- 


40 


- 20 


^ 


' ' Gh'prediction ' ' 1 ' - 
^ — Scintillation / - 




SO 


110 
- 


- 


\. , , J i 




: 




, I , , ' t^S-r-T-l , , r 





0.9 

i . 


0.92 0.94 0.96 0.98 

, 1 , . , 1 


V 


1 ~ 



0.4 0.6 

Sin(inclination) 



FIG. 2: Current constraints on the orbital inclination angle (i) 
of the PSR J1141— 6545 system. The smooth curve represents 
values of A^ 2 obtained when sin i was held fixed at a range of 
values from 0.2 to 1 (i.e. i from 10° to 90°) and the companion 
mass at the value determined by mtot = 2.2892 Mq and the 
mass function (f p ) = 0.176701 M . 

i = 76° ± 2°. 5 derived from the orbital modulation of the 
pulsar's scintillation velocity 13]. 

Figure[2]shows a goodness-of-fit plot for sin i values be- 
tween 0.2 and 1, in which the companion mass was set to 
a value determined by the precisely known systemic mass 
mtot and the pulsar mass function, f p = (m c sim) 3 /mj? ot 
= 0.176701 M . Such an approach is justifiable as the 
relativistic to is determined at a very high precision and 
hence traces out a unique locus in the mass-mass dia- 
gram, with each point on the locus implying a unique 
inclination angle. While the use of the DDGR-derived 
mtot means assumption of GR, the timing model used is 
DD. The two-sigma range in the inclination angle (72°. 6 
to 79°. 4) indicated by this curve is consistent with com- 
pletely independent constraints derived from the scintil- 
lation measurements and thus, effectively, provides yet 
another confirmation of GR in this system. 

The revised mass estimate of PSR J1141— 6545 makes 
it one of the lighter neutron stars currently known. Al- 
though the pulsar's progenitor was initially the lighter 
star in the original binary pair, mass transfer from the 
white dwarf's progenitor made it large enough to explode 
as a supernova even though its companion was not large 
enough to do so [lj|. Being the second star to evolve 
in the binary, PSR J1141— 6545 is somewhat similar to 
pulsar B in the double pulsar system and the companion 
to PSR J1756-2251 %^. It seems probable therefore, 
that the progenitors to these pulsars were of lower mass 
and that the mass of a neutron star is related to that of 
the progenitor. Recycled pulsars, or those with higher 
eccentricities, appear to be heavier [2fJ, [21j either due to 
mass accretion [22] or higher initial masses. 

Alternative Theories of Gravity. — Deviations from GR 
inherent to most alternative theories of gravity are 
best detectable in strong gravitational fields, and hence 
binary-pulsar measurements are indispensable in test- 
ing such theories [23|]. In particular, the gravitation- 
ally asymmetric nature of the PSR J1141— 6545 system 
makes it a unique testing ground for tensor-scalar theo- 



ries Q, the best motivated alternatives to GR. This is 
because for such systems, tensor-scalar theories predict 
the emission of a large amount of dipolar scalar waves 
in the form of gravitational radiation loss (as opposed to 
predominantly quadrupolar radiation predicted by GR). 
This effect is strongly suppressed in double neutron star 
binaries, and as a result our timing measurements of 
PSR J1141— 6545 impose much tighter constraints on 
these theories than those possible by PSRs B1913+16, 
B1534+12 and J0737-3039 [IGjJ. In fact, the measure- 
ment of an orbital decay that is in accord with GR would 
naturally exclude a significant part of the theory space. 
A parameter of particular interest is the linear coupling 
strength ao , for which the asymptotic limit scales as the 
square root of the measurement uncertainty in P 8]. 
This limit corresponds to strongly non-linear coupling 
(the quadrature coupling strength, (3q ~10 or larger). A 
simple extrapolation based on our current measurements 
yields a^ < 3.4 x 1CP 6 in this regime. This is three times 
smaller than the Cassini bound (a^ < 1-15 x 10" 5 ), and 
thus provides the current best limit in that part of the 
theory space. For smaller values of (3q (i.e. weakly non- 
linear coupling) however, the limit is weaker and is a func- 
tion of /3o ■ In the special case of linear coupling (Brans- 
Dicke; /3o = 0), it is reduced by a factor of 4er^ s , where 
e ns is the self-gravity of the neutron star Q ■ Assuming 
e ns = 0.2, we derive a§ < 2.1 x 10~ 5 for this regime, which 
is nearly twice the Cassini limit. Thus, the limit imposed 
by our current data is weaker than that from Cassini 
for linear or weakly non-linear coupling but becomes in- 
creasingly stringent for strongly non-linear coupling and 
it should improve upon the Cassini bound somewhere 
within the range 5 £> (3q S 10. As our timing precision im- 
proves over time, these limits are expected to become 
even more stringent, eventually cutting well below the 
Cassini bound at much lower values of /3o ■ 

Conclusions and Future Prospects. — Our long timing 
campaign on PSR J1141— 6545 has led to precise mea- 
surements of its PK parameters d) , 7, and P\, , providing 
strong confirmation of GR in gravitationally asymmet- 



ric binary systems. The measured orbital decay is in 
agreement with the GR prediction and the Shapiro delay 
is consistent with independent measurements of the or- 
bital inclination angle. The pulsar and companion masses 
are determined to be m p = 1.27 ± 0.01 M© and m c — 
1.02 ± 0.01 M Q respectively. Stringent limits are placed 
on tensor-scalar theories of gravity, with an asymptotic 
limit of ojq < 3.4 x 10 -6 , cutting well below the Cassini 
limit for large values of the coupling parameter /3 and 
hence providing the lowest limit in that range. For small 
values of /?o , the limit is a^ < 2.1 x 10~ 5 and thus weaker 
than that derived from the Cassini experiment. 

Continued timing of this pulsar looks very promising. 
Currently Pa is determined at a 6% precision, the uncer- 
tainty in which is expected to decrease as T~ 5 / 2 , where 
T is the observing time span. We thus anticipate ap- 
proaching a precision near 2% by 2012. At such high pre- 
cision, contaminations from the kinematic and Galactic 
contributions will start dominating the error budget [18j • 
While the Galactic term P^Jai j g narc j ^o determine accu- 
rately, the kinematic Shklovskii contribution P^ ln can be 
assessed better when an independent measurement be- 
comes available for the pulsar's transverse motion. Fortu- 
nately, the pulsar's location is such that P^» al (dominated 
by differential rotation in the plane of the Galaxy) is of 
opposite sign to P£ m , and thus may cancel out a large 
fraction of it if the pulsar is moving at ~100 kms _1 or 
faster, potentially resulting in a net contamination that 
is well below 1%. This will enable even more precise tests 
in the future and will place the most stringent constraints 
on alternative theories of gravity. 
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